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We consider a dynamical system of phantom scalar field under exponential potential in background 
of loop quantum cosmology. In our analysis, there is neither stable node nor repeller unstable node 
but only two saddle points, hence no Big Rip singularity. Physical solutions always possess potential 
energy greater than magnitude of the negative kinetic energy. We found that the universe bounces 
after accelerating even in the domination of the phantom field. After bouncing, the universe finally 
enters oscillatory regime. 
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INTRODUCTION 



Recently, present accelerating expansion of the uni- 
verse has been confirmed with observations via cosmic 
microwave background anisotropics [J, large scale 
galaxy surveys [3] and type la supernovae 0, H[ . How- 
ever, the problem is that the acceleration can not be un- 
derstood in standard cosmology. This motivates many 
groups of cosmologists to find out the answers. Propos- 
als to explain this acceleration made till today could be, 
in general, categorized into three ways of approach @. 
In the first approach, in order to achieve acceleration, 
we need some form of scalar fluid so called dark energy 
with equation of state p = wp where w < —1/3. Various 
types of model in this category have been proposed and 
classified (for a recent review see Ref. The other 

two ways are that accelerating expansion is an effect of 
backreaction of cosmological perturbations [9[ or late ac- 
celeration is an effect of modification in action of gen- 
eral relativity. This modified gravity approach includes 
braneworld models (for review, see [10]). Till today there 
has not yet been true satisfied explanation of the present 
acceleration expansion. 

Considering dark energy models, a previous first year 
WMAP data analysis combined with 2dF galaxy sur- 
vey and SN-Ia data and even a previous SN-Ia analy- 
sis alone favor w < — 1 than cosmological constant or 
quintessence [ll|, [l2| . Precise observational data analy- 
sis (combining CMB, Hubble Space Telescope, type la 
Supernovae and 2dF datasets) allows equation of state 
p = wp with constant w value between -1.38 and -0.82 
at the 95 % of confident level [13]. The recent WMAP 
three year results combined with Supernova Legacy Sur- 
vey (SNLS) data when assuming flat universe yields 
— 1.06 < w < —0.90. However without assumption 
of flat universe but only combined WMAP, large scale 



structure and su per nova data implies strong constraint, 
1141 ]. While assuming flat universe, the 
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first result from ESSENCE Supernova Survey la com- 
bined with SuperNova Legacy Survey la gives a con- 
straint oiw = — 1.07 ± 0.09 [la]. Interpretation of various 
data brings about a possibility that dark energy could be 
in a form of phantom field-a fluid with w < — 1 (which 
violates dominant energy condition, p > \p\) rather than 
quintessence field [la, UJ|, LL8[ . The phantom equation of 
state p < —p can be attained by negative kinetic energy 
term of the phantom field. However there are some types 
of braneworld model pjl as well as Brans-Dicke scalar- 
tensor theory [20[ and gravitational theory with higher 
derivatives of scalar field [21] that can also yield phan- 
tom energy. There has been investigation on dynamical 
properties of the phantom field in the standard FRW 
background with exponential and inverse-power law po- 
tentials byJ^Jg^Jp, ^ and with other forms of po- 
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tential by [25|, [26|, [27| . These studies describe fates of the 
phantom dominated universe with different steepness of 
the potentials. 

A problem for phantom field dark energy in standard 
FRW cosmology is that it leads to singularity. Fluid with 
w less than -1 can end up with future singularity so called 
the Big Rip [28[ which is of type I singularity according 
to classification by [29|, |30[ . The Big Rip singularity cor- 
responds to a — > oo, p — > oo and \p\ — > oo at finite time 
t — > t s in future. Choosing particular class of potential 
for phantom field enables us to avoid future singularity. 
However, the avoidance does not cover general classes of 
potential [2f|. In addition, alternative model, in which 
two scalar fields appear with inverse power-law and ex- 
ponential potentials, can as well avoid the Big Rip singu- 
larity [31]. The higher-order string curvature correction 
terms can also show possibility that the Big-Rip singu- 
larity can be absent [32[. 

Since phantom dominated FRW universe possesses sin- 
gularity problem as stated above, in this work, instead of 
using standard FRW cosmology, the fundamental back- 
ground theory in which we are interested is Loop Quan- 
tum Gravity-LQG. This theory is a non-perturbative 
type of quantization of gravity and is background- 
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independent [33|, [34| . It has been applied in cosmological 
context as seen in various literatures where it is known 
as Loop Quantum Cosmology-LQC (for review, see Ref. 
[35l]). Effective loop quantum modifies standard Fried- 
mann equation by adding a correction term —p 2 /p\ c into 
the Friedmann equation [H, Ezl EH, [H, Ell- When this 
term becomes dominant, the universe begins to bounce 
and then expands backwards. LQG can resolve of singu- 
larity problem in various situations [H, [3?], EH EH • How- 
ever, derivation of the modified term is under a condition 
that there is no matter potential otherwise, in presence of 
a potential, quantum correction would be more compli- 
cated [43]. Nice feature of LQC is avoidance of the future 
singularity from the correction quadratic term —p 2 /p\ c in 
the modified LQC Friedmann equation [44] as well as the 
singularity avoidance at semi-classical regime [45]. The 
early-universe inflation has also been studied in the con- 
text of LQC at semi-classical limit 0, El, Ell EE HI HI • 
We aim to investigate dynamics of the phantom field and 
its late time behavior in the loop quantum cosmological 
context, and to check if the loop quantum effect could 
remove Big Rip singularity from the phantom dominated 
universe. The study could also reveal some other inter- 
esting features of the model. 

We organize this article as follows: in section [III we in- 
troduce LQC Friedmann equation, after that we briefly 
present relevant features of the phantom scalar field in 
section IIIII Section IIVI contains dynamical analysis of 
the phantom field in LQC background with exponential 
potential. The potential is a simplest case due to con- 
stancy of its steepness variable A. Two real fixed points 
are found in this section. Stability analysis yields that 
both fixed points are saddle points. Numerical results 
and analysis of solutions can be seen in section IVl where 
we give conditions for physical solutions. Finally, conclu- 
sion is in section [VTl 



where c and p are respectively conjugate connection and 
triad satisfying {c,p} = j/3Mp. Dot symbol denotes 
time derivative. These are two variables in the simpli- 
fied phase space structure under FRW symmetries [351 ] . 
Here Mp = (87rG) _1 is square of reduced Planck mass, 
G is Newton's gravitational constant and 7 is Barbero- 
Immirzi dimensionless parameter. There are relations 
between the two variables to scale factor as p = a 2 and 
c = 7a. The parameter p is inferred as kinematical length 
of the square loop since its order of magnitude is simi- 
lar to that of length. The area of the loop is given by 
minimum eigenvalue of LQG area operator. C m is the 
corresponding matter Hamiltonian. Using the Eq. ([2]) 
with constraint from realization that loop quantum cor- 
rection of effective Hamiltonian C e ff is small at large scale, 
C e ff ~ [H, [H, Hill, E3 , one can obtain (effective) modi- 
fied Friedmann equation in flat universe: 
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(3) 



where p\ c — \f3 / (Wttj 3 G 2 h) is critical loop quantum 
density, h is Planck constant and p t is total density. 



III. PHANTOM SCALAR FIELD 

The energy density p and the pressure p of the phantom 
field contain negative kinetic term. They are given as [l6[ 



The conservation law is 



p + 3H(p + p) =0. 



(4) 
(5) 

(6) 



II. LOOP QUANTUM COSMOLOGY 

LQC naturally gives rise to inflationary phase of the 
early universe with graceful exit, however the same mech- 
anism leads to a prediction that present-day acceleration 
must be very small [46[. At late time and at large scale, 
the semi-classical approximation in LQC formalisms can 
be validly used [5l|. The effective Friedmann equation 
can be obtained by using an effective Hamiltonian with 
loop quantum modifications [38|, [44], [52[ : 



eff 



3M 2 
7 2 /2 2 



a sin 2 (/ic) + C n 



(i) 



The effective constraint (pQ) is valid for isotropic model 
and if there is scalar field, the field must be free, mass- 
less scalar field. The equation (pQ), when including field 
potential, must have some additional correction terms 
[431 ]. In this scenario, the Hamilton's equation is 



P = {p,C eff } = 



7 dCeff 
~3M 2 dc 



(2) 



Using the Eqs. dH), J5]) and ©, we obtain Klein-Gordon 
equation: 



3iJ0 - V = , 



(7) 



where V = d V/ dcf) and the negative sign comes from the 
negative kinetic terms. The phantom equation of state 
is therefore given by 



w = 



V 



2V 



2V 



(8) 



From the Eq. ([8]), when the field is slowly rolling, as long 
as the approximation, <p 2 ~ holds, the approximated 
value of w is -1. When the bound, (j) 2 < 2V holds, w is 
always less than -1. 

As mentioned before in sections [J and [III there has not 
yet been a derivation of effective LQC Friedmann equa- 
tion in consistence with a presence of potential. Even 
though, the Friedmann background is valid only in ab- 
sence of field potential, however, investigation of a phan- 
tom field evolving under a potential is a challenged task. 
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Here we also neglect loop quantum correction effect in 
the classical expression of Eqs. dU) and ([5]) (see Refs. 
(43^ and [5^ for discussion). 



IV. DYNAMICAL ANALYSIS 

Differentiating the Eq. (j3]) and using the fluid Eq. (j6]), 
we obtain 



fT _ (P + P) (, 2p 
2M 2 \ Pic 



(9) 



The Eqs. (|3]) , (|6]) and ([9]), in domination of the phantom 
field, become 



p = -3Hp 1 



t> 2 + 2V \ 
h 2 -2V ' 



2M 2 V* Pk 



(11) 



(12) 



We define dimensionless variables following the style of 

M 

* = , Y=-^-, Z=±, (13) 

VQM p H V3M P H pic 

M P V „ V V" d Id 
^ ^> 2, ^^^37, (14) 



1/ 



(V) dN Hdt 



where N = In a is e-folding number. Using new variables 
in Eqs. (jHJ) and (fTUj) . the equation of state is rewritten 



as 



X 2 + Y 2 
X 2 - Y 2 ' 



(15) 



where |X| 7^ |Y"| and the Friedmann constraint is reex- 
pressed as 



(-x 2 + y 2 )(i-z) = i. 



(16) 



Clearly, if \X\ ^ \Y\, following the Eq. (JTSJ), then Z ^ 1. 
Using the new defined variables above, Eq. ([T2|) becomes 



3X 2 (1 - 2Z) . 



The acceleration condition, 



- = H + H z > 0, 
a 



(17) 



(18) 



in expression of the new variables, is therefore 

3X 2 (2Z-1)<1. (19) 



Divided by the Eq. ([16]) , the acceleration condition under 
the constraint is 



1-2Z 



1 - (Y 2 /X 2 ) \1-Z 



<1, 



(20) 



where the conditions |X| 7^ \Y\ and Z / 1 must hold. 
As we consider Z = p/p\ c with p = — (</> 2 /2) + V, we can 
write 



3M 2 # 2 



-x 2 + y 2 . 



(21) 



With the condition |X| 7^ |Y|, clearly from Eq. ([2T]h we 
have one additional condition, 



A. Autonomous system 



Differential equations in autonomous system are 

^ = -3X - J\ XY 2 - 3X 3 (1 - 2Z) , (22) 



dr 

diV 
dZ 
dN 
dX 
dN 



XXY - 3X 2 Y (1 - 2Z) , (23) 



-3Z 1 



x 2 + r s 
x 2 - r 2 y ■ 

-^(r-i)A 2 x. 



Here we will apply exponential potential, 
V{4>)=V e W {—±-4>), 



(24) 
(25) 

(26) 



to this system. The potential is known to yield power- 
law inflation in standard cosmology with canonical scalar 
field. Its slow-roll parameters are related as e = rj/2 = 
1/P where A = j^fP and P > 1 [5j| E^. Although 
the potential has been applied to the quintessence scalar 
field with tracking behavior in standard cosmology [57}, 
the quintessence field can not dominate the universe due 
to constancy of the ratio between densities of matter and 
quintessence field (see discussion in Ref. [7]). In case 
of phantom field in standard cosmology under this po- 
tential, a stable node is a scalar-field dominated solution 
with the equation of state, w = -1 - A 2 /3 0, [III, EH . In 
our LQC phantom domination context, from Eq. ([25]) . 
we can see that for the exponential potential, T = 1. 
This yields trivial value of dX/dN and therefore A is a 
non-zero constant otherwise the potential is flat. 



1 The relation = p/SH 2 M 2 = -X 2 + Y 2 = 1 can not be 
applied here since it is valid only for standard cosmology with 
flat geometry. 



B. Fixed points 

Let / = dX/dN,g = dY/dN and h = dZ/dN. We 
can find fixed points of the autonomous system under 
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FIG. 1: Three-dimensional phase space of X,Y and Z. The 
saddle points (a) (-0.40825, 1.0801, 0) and (b) (-0.40825, - 
1.0801, 0) appear in the figure. A is set to 1. In region Z < 0, 
the solutions (red and blue lines) are non physical. In this 
region, Z — > — oo when (X, Y) — > (0, 0). The green lines (class 
I) are in region \X\ > \Y\ and Z > 1 but they are also non 
physical since they correspond to imaginary H values. The 
only set of physical solutions (class II) is presented with black 
lines. They are in region \Y\ > \X\ and range < Z < 1. 
This is the region above (a) and (b) of which H takes real 
value. There are separatices \X\ = \Y\, Z = and Z = 1 in 
the system (see section |VB]) . 



condition: 

(/, 9, h) |(x c> y c ,z c )=0. 
The are two real fixed points of this system: 2 



• Point (a) : (-^, 



A 2 

1 + T ,0), 



• Point (b) : + 0). 



V6' 



(27) 

(28) 
(29) 



C. Stability Analysis 

Suppose that there is a small perturbation SX, SY 
and 8Z about the fixed point (X C ,Y C , Z c ), i.e., 

X = X C + SX, Y = Y C + SY, Z = Z C + 5Z. (30) 



2 The other two imaginary fixed points (i,0, 0) and (—2,0,0) also 
exist. However they are not interesting here since we do not 
consider model that includes complex scalar field. 



FIG. 2: Phase space of the kinetic part X and potential part 
Y (top view). The saddle points (a) (-0.40825, 1.0801) and 
(b) (-0.40825, -1.0801) are shown here. The blue lines and 
red lines are in the region Z < which is non physical. Green 
lines are of class I solutions which yields imaginary H. Only 
class II solutions shown as black lines are physical with real 
H value. 



From Eqs. ([22]), ([23} and ([24j), neglecting higher order 
terms in the perturbations, we obtain first-order differ- 
ential equations: 



dN 



SX \ 


SX 


SY 


= M\ SY 


SZ 


\ SZ 



(31) 



The matrix Ai defined at a fixed point (X C ,Y C , Z c ) is 
given by 

Af= I I f (32) 

I M K M / 

V dX dY dZ / (X=X C ,Y=Y C ,Z=Z C ) 

We find eigenvalues of the matrix M. for each fixed point: 
• At point (a): 

« A 2 



• At point (b): 



jii = X 



112 



-A 



-A 



Ms 



A 2 
2 



(33) 



(34) 



From the above analysis, each point possesses eigenval- 
ues with opposite signs, therefore both point (a) and (b) 
are saddle. Results from our analysis are concluded in 
TABLE [U Location of the points depends only on A and 
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Name 


X 


Y 


Z 


Existence 


Stability 


w 


Acceleration 


(a) 


X 

V6 







All A 


Saddle point for all A 


_i_ a! 

X 3 


For all A (i.e. A 2 > -2) 


(b) 


A 

V6 







All A 


Saddle point for all A 


\ 2 

-1 — — 

1 3 


For all A (i.e. A 2 > -2) 



TABLE I: Properties of fixed points of phantom field dynamics in LQC background under the exponential potential. 




-1.5 -1 -0.5 0.5 1 1.5 



X 



FIG. 3: Phase space of the kinetic part X and potential part 
Y in standard general relativist ic case. The location of points 
(a) and (b) in Fig. [2] are on the trajectory solutions here. 
This plot shows dynamics of phantom field in standard cos- 
mological background without any other fluids. In presence 
of a barotropic fluid with any equation of state, the point (a) 
and (b) correspond to the Big Rip [23l [25j. 



the points exist for all values of A. Both points corre- 
spond to the equation of state —1 — A 2 /3, that is to say, 
it has phantom equation of state for all values of A ^ 0. 
Since there is no any attractor in the system, a phase tra- 
jectory is very sensitive to initial conditions given to the 
system. The stable node (the Big Rip) of the standard 
general relativistic case in presence of phantom field and 
a barotropic fluid, disappears here (see @). 



V. NUMERICAL RESULTS 

Numerical results from the autonomous set ([22]) , ([23]) 
and ([24]) are presented in Figs. [1] and [2] where we set 
A = 1. Locations of the two saddle points are: point 
(a) (X c = -0.40825, Y c = 1.0801, Z c = 0) and point (b) 
(X c = -0.40825, Y c = -1.0801, Z c = 0) which match 
our analytical results. In Fig. O we present a trajectory 
solution of a phantom field evolving in standard cosmo- 
logical background for comparing with the trajectories in 
Fig. [2] when including loop quantum effects. The stan- 



dard case has only simple two trajectories corresponding 
to a constraint — X 2 -\-Y 2 — 1. This is attained when tak- 
ing classical limit, Z = 0. In loop quantum case, since 
there is no any stable node and the solutions are sen- 
sitive to initial conditions, we need to classify solutions 
according to each domain region separated by separatri- 
ces \X\ = |y|, Z = and Z = 1 so that we can analyze 
them separately. Note that the condition, Z > must 
hold for physical solutions since the density can not be 
negative or zero, i.e. p > 0. The blue lines and red lines 
in Figs. [T]and [2] are solutions in the region Z < hence 
are not physical and will no longer be of our interest. 
From now on we consider only the region Z > 0. In 
regions with \X\ > \Y\, the solutions therein are green 
lines (hereafter classified as class I). The other regions 
with \Y\ > \X\ contain solutions seen as black line (clas- 
sified as class II). Note that all solutions can not cross 
the separatices due to conditions in Eqs. ([T6]h ([20]) and 

A. Class I solutions 

Consider the Friedmann equation ([TO]) , the Hubble pa- 
rameter, H takes real value only if 

^(-H^)*"- (35) 

Divided by H 2 on both sides, the expression above be- 
comes 

(-X 2 + F 2 )(l-Z) > 0. (36) 

It is clear from ([36]) that, in order to obtain real value 
of H, class I solutions (green line) must obey both con- 
ditions \X\ > \Y\ and Z > 1 together. However, when 
imposing \X\ > \Y\ to the Eq. ([21]), we obtain Z < 
instead. This contradicts to the required range Z > 1. 
Therefore this class of solutions does not possess any real 
values of H and hence not physical solutions. 

B. Class II solutions 

Proceeding the same analysis done for class I, we found 
that in order for H to be real, class II solutions must 
obey both \Y\ > \X\ and < Z < 1 together. Moreover 
when imposing \Y\ > \X\ into Eq. ([2T]h we obtain Z > 0. 
Therefore as we combine both results, it can be concluded 
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FIG. 4: Evolution of H with time of a class II solution. Set 
values are A = l,pic — 1.5, Vb — 1 and Mp — 2. The uni- 
verse undergoes acceleration from the beginning until reach- 
ing turning point at p — pi c /2 — 0.75 where H — if max — 
0.17678. Beyond this point, the universe expands with decel- 
eration until halting (H = 0) at p w p\ c = 1.5. After halting, 
it undergoes contraction until H bounces. The oscillating in 
H goes on forever. 



that class II solutions can possess real H value in the 
region \Y\ > \X\ and < Z < 1, i.e. < p < p\ c . 
The bound is slightly different from the case of canonical 
scalar field in LQC (see Ref. of which the bound is 
< p < pic- The class II is therefore the only class of 
physical solutions. 

For class II solutions, we consider another set of au- 
tonomous equations from which the evolution of cosmo- 
logical variables are conveniently obtained by using nu- 
merical approach. In the new autonomous set, instead 
of using A 7 ", which could decrease after the bounce from 
LQC effect, time is taken as independent variable. We 
define new variable as 

= 5. (37) 
The Eqs. and ([T2]) are therefore rewritten as 



H 

S 



s 2 



2M 2 L 
-3HS 



2 

Pic 



2 



y(0) 



(38) 
(39) 



The Eqs. (|37j), (J38J) and ([39]) form another closed au- 
tonomous system. Numerical integrations from the new 
system yield result plotted in Figs. [4] and [5] in which set 
values are A = l,pi c = 1.5, Vb = 1 and Mp = 2. From 
Eq. Q the slope of H with respect to p, diJ/dp, is fiat 
when p = pic/2 |59| . Another fact is 



'd 2 H 
K dp 2 



p=p\ c /2 



Mp y/3p 



< o, 



(40) 



FIG. 5: Time evolution of potential energy density (P.E.), 
kinetic energy density (K.E.) and p — K.E. + P.E. of the 
field for a class II solution. K.E. is always negative and, at 
late time, it goes to — oo while P.E. is always positive, p is 
maximum when p w p\ c — 1.5. Other features are discussed 
as in Fig. |4] 




FIG. 6: Oscillation in kinetic energy density (K.E.) that con- 
tributes to oscillation in p. This is a zoom-in portion of the 
Fig. [5] 



hence, as p = p\ c /2, H takes maximum value, i^ max = 
yj p\ c /12Mp. This result is valid in LQC scenario regard- 
less of types of fluid. In Figs. [H and [5l with set param- 
eters given above, as p = p\ c /2 = 0.75, H is maximum, 
F max = 0.17678. When H « 0, i.e. p is approximately 
pic = 1.5, the expansion halts and then bounces. At 
this bouncing point, the dynamics enters loop quantum 
regime which is a quantum gravity limit. Beyond the 
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bounce, H turns negative, i.e. contracting of scale fac- 
tor. The universe undergoes accelerating contraction un- 
til reaching H m { n . After that it contracts deceleratingly 
until bouncing at H w 0. The universe goes on faster 
bouncing forward and backward. The faster bounce in 
H is an effect from the faster bounce in p as illustrated in 
Fig. [5] where the red line represents potential energy den- 
sity V ((/)), the black line represents kinetic energy density 
— (j) 2 /2 and the blue line is total energy density p. Oscilla- 
tion in p is from oscillation in the field speed <j) and there- 
fore oscillation in K.E. as shown in Fig. [6l This hence 
contributes to oscillation in p. The negative magnitude 
of kinetic energy density becomes larger and larger as the 
field rolling faster and faster up the potential. The expo- 
nential potential energy density therefore becomes larger 
and larger. This results in oscillation of p and affects in 
oscillation of H about the bounce H = 0. With a differ- 
ent approach, recently a similar result in H oscillation is 
also obtained by Naskar and Ward [60[ . 

VI. CONCLUSION 

A dynamical system of phantom canonical scalar 
field evolving in background of loop quantum cos- 
mology is considered and analyzed in this work. 
Exponential potential is used in this system. Dy- 
namical analysis of autonomous system renders 
two real fixed points (— A/\/6, y^l + A 2 /6 , 0) and 

(-A/V / 6,- v /l + A2 /6, 0), both of which are saddle 
points corresponding to equation of state, w = — 1 — A 2 /3. 
Note that in case of st andard co smology, the fixed point 
(X c , Y c ) = (-A/V6, v /l + A 7 6 ) is the Bi S Ri P attractor 
with the same equation of state, w = — 1 — A 2 /3 [24^ . 
Due to absence of stable node, the late time behavior 
depends on initial conditions given. Therefore we do 
numerical plots to investigate solutions of the system 
and then classify the solutions. Separatrix conditions 
\X\ |F|, Z 1 and Z / arise from equation of 
state ([15]) . Friedmann constraint ([T6]) and definition of 
Z in Eq. ([2T]) . At first, we consider solutions in region 



Z > 0, i.e. p > for physical solutions. Secondly, 
within this Z > region, we classify solutions into class 
I & II. Solutions in region \X\ > \Y\ and Z > 1 are of 
class I. However, in order to obtain real value of H in 
class I, Z must be negative which contradicts to Z > 1. 
Therefore the class I solutions are non physical. Class II 
set is identified by \Y\ > \X\ and < Z < 1. It is an 
only set of physical solutions since it yields real value of 
H. In class II set, the universe undergoes accelerating 
expansion from the beginning until p = p\ c /2 where 
H = i^ m ax = \J pic/12Mp. After that the universe 
expands deceleratingly until it bounces, i.e. stops 
expansion H « at p « p\ c . At the bounce the universe 
enters quantum gravity regime. Contraction with 
backward acceleration happens right after the bounce, 
however the contraction does not go on forever. When 
the universe reaches minimum value of negative H, the 
contraction turns decelerated, i.e. contracts slower and 
slower down. The universe, after undergoing contraction 
to minimum spatial size, bounces again and starts to 
expand acceleratingly. Our numerical results yield that 
oscillation in H becomes faster as time passes. 
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